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Abstract 



We consider a late-time cosmological model based on a recent proposal that 
the infinite-bare-coupling limit of superstring/M-theory exists and has good 
phenomenological properties, including a vanishing cosmological constant, and 
a massless, decoupled dilaton. As it runs away to +oo, the dilaton can play 
the role of the quintessence field recently advocated to drive the late-time 
accelerated expansion of the Universe. If, as suggested by some string theory 
examples, appreciable deviations from General Relativity persist even today 
in the dark matter sector, the Universe may smoothly evolve from an ini- 
tial "focusing" stage, lasting untill radiation-matter equality, to a "dragging" 
regime, which eventually gives rise to an accelerated expansion with frozen 
fi(dark energy) /O (dark matter). 
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I. INTRODUCTION 



According to recent astrophysical observations, our Universe, since a red shift of (9(1), ap- 
pears to have undergone a phase of accelerated expansion |j],[2|] . This result can be combined 
with the recent estimates of the average mass density of the Universe ||, Q m ~ 0.3 — 0.4 
(in critical units), and with recent measurements of the CMB anisotropy peaks ||, pointing 
at a nearly critical total energy density, — 1- One is then led to the conclusion that 
the present cosmological evolution, when described in terms of an effective fluid entering 
Einstein's equations, should be (marginally) dominated by a "dark energy" component p x 
characterized by a (sufficiently) negative effective pressure, p x < —p x /3. 

The simplest candidate for such a missing energy is a positive cosmological constant A, of 
order Hq. Such an identification, however, unavoidably raises a series of difficult questions. 
In particular: a) Why is A so small in particle physics units? Explaining a finite but very 
small value for A may turn out to be even harder than finding a reason why it is exactly zero. 
This is the so-called fine-tuning problem for A, see for instance [[|; and b) Why is A ~ p m0 , 
where p m o is the present value (in Planck units) of the (dark) matter energy density? This 
is the so-called "cosmic coincidence" problem |J . 

At present, the most promising scenarios for solving (at least part of) the above prob- 
lems introduce a single scalar field, dubbed "quintessence" J7j , whose potential goes to zero 
asymptotically (leaving therefore just the usual puzzle of why the "true" cosmological con- 
stant vanishes). The scalar field slowly rolls down such a potential reaching infinity (and 
zero energy) only after an infinite (or very long) time. While doing so, quintessence pro- 
duces an effective, time-dependent, cosmic energy density p x accompanied by a sufficiently 
negative pressure, i.e. a sort of effective cosmological constant. By making A e ff ~ H 2 time- 
dependent, this can naturally explain the smallness of the present effective vacuum energy 
density. However, if, as in General Relativity, dust energy and an effective cosmological 
constant have different time dependence, it can hardly explain why A ~ p m0 . For a recent 
review of the relative merits of a cosmological constant and quintessence, see Ref. [§]. 

As far as identifying quintessence is concerned, the inflaton itself could be a candidate 
||. But also more exotic possibilities have been considered, in particular some motivated by 
the wish to solve the above-mentioned cosmic coincidence problem fT0] , [TT|j . In any case, as is 
the case for the inflaton, the quintessence field does not have, as yet, an obvious place in any 
fundamental theory of elementary particles. One should also mention, at this point, that, if 
quintessence may help with the problems typical of the cosmological constant interpretation, 
it is likely to create a new one of its own: in order to play its role, the quintessence field must 
be extremely light and can thus mediate a new long-range (of order Hq 1 ) force, which is 
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strongly constrained observationally. This is an important constraint to be imposed on any 
specific scalar field model of quintessence, either minimally or non-minimally [T2|] coupled to 
gravity. 

At first sight, the search for a quintessence candidate in particle physics looks easier 
than the one for an inflaton. For instance, fundamental or effective scalar fields with po- 
tentials running to zero at infinity are ubiquitous in supersymmetric field theories and/or in 
string/M-theory. They are usually referred to as moduli fields since, in perturbation theory, 
they parametrize the space of inequivalent vacua and correspond to exactly flat directions 
(equivalently, to exactly massless fields). Non-perturbative effects (e.g. gauge-theory instan- 
tons) are expected to lift these flat directions, just preserving those that correspond to small 
or vanishing coupling. Examples are the run-away vacua of supersymmetric gauge theories 
(see, for instance, |13| for quintessence models based on the latter possibility), or the dilaton 



modulus (p in the limit — > — oo. 

However, if we were to take one of these moduli as quintessence, we would immediately 
run into the problem that the acceleration of the Universe should be accompanied by a drift 
of interactions towards triviality. For Newton's constant, and even more so for the fine- 
structure constant, this kind of time variations is very strongly constrained. Furthermore, 
typical couplings of moduli fields to ordinary matter are of gravitational order, and this 
creates the already mentioned problem of new, unwanted long-range forces. For all these 



reasons the conventional attitude towards moduli fields has been (see e.g. [|14|) to postulate 
that they develop non-perturbative potentials, providing them with both a mass and a 
freezing mechanism (see however |H| for an alternative that is closer, in spirit, to the one 
advocated here). 

Another possible problem with the identification of a string modulus with quintessence 
is that we would like to freeze the moduli at values that provide the correct values of the 
coupling constant and unification scale of grand unified theories (GUTs). For instance, the 



dilaton and compactification volume Vq should be frozen at values such that |T6] 

a- G l UT ~ (M P /M S ) 2 ~ eT* , M P /M GUT ~ a G %l g x J\ (1.1) 

where Mp, M s and Mqut are the Planck, string and GUT scales, g s = e^V^M® is the string 
coupling, and is the tree-level effective four-dimensional coupling (thus, in more standard 



string-theory notation [|Tj, our dilaton is related to the real part of the S modulus by Rej,?} 

Unfortunately, it looks unlikely that non-perturbative effects will be significant enough 
in this region to stabilize the moduli. Also, perturbative unification gives too low a value for 
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M P /M GUT [|T^,|T^]. In this respect, the situation can be drastically improved by considering 
the M-theory limit, g s — > oo, while still keeping the four-dimensional effective couplings 
perturbative (5> 1) flSfl . Even then, the moduli would presumably freeze out in a typical 
(and cosmologically tiny) particle-physics time, and therefore cannot implement the conven- 
tional, slow-roll quintessential scenario. In spite of these difficulties, unconventional models 
of quintessence based on the stabilization of the dilaton in the perturbative regime are not 



completely excluded, as recently discussed by one of us 

There is, however, another possibility for making the dilaton a candidate for quintessence. 
As we have already mentioned, the region of large negative corresponds to the trivial 
vacuum. The idea that the Universe may have started, long before the big bang, in this 
region is actually the basis of the so-called pre-big bang scenario in string cosmology (for 



recent reviews see 0]). Here we are asking instead whether the dilaton can play the role of 
quintessence at very late times (such as today), not by evolving towards — oo and triviality, 
but by going towards +oo and strong coupling. Such a proposal looks absurd at first since, 
if we do not see a drift towards zero coupling, we do not experience one towards increasing 
strength either. In order for this idea to make sense we have to assume that the strong- 
coupling limit of string/M-theory exists, is smooth, and resembles our world. Can this make 
sense at all? 



It has been argued by one of us |?TJ that the answer to the above question can be 
affirmative, if we assume a certain structure of the quantum loop corrections to the string 
effective action suggested by large- N counting arguments. In the strong coupling limit (which 
could either be the self-dual value S — 1 or, if ^-duality is broken, (f> — > +oo) gravitational 
and gauge coupling would be determined entirely by loop corrections (as in the so-called 
induced-gravity /gauge idea f22|), and would "saturate" at "small" values because of the 
large number of fields entering the loops (e.g. the large number N of gauge bosons, or the 
large value of the quadratic Casimir Ca, for gauge groups like E$). Typically, eq. (|TTT| ) 
would be replaced (at <fi ^> 1) by 

a- G 1 UT ~C A + 0(e-' t> ) , (M P /M s ) 2 ~N + 0(e-*). (1.2) 

In this picture there is naturally an asymptotic decoupling mechanism of ordinary matter 
to the dilaton, whose effective mass goes to zero at late times. The problem remains, of 
course, of explaining why the cosmological constant vanishes in superstring/M-theory, not 
only at zero coupling where supersymmetry protects it, but also at infinite (bare) coupling. 
Possibly, some new, stringy symmetry can explain this. It will simply be assumed to be the 
case in this paper. 

As the dilaton is non-universally coupled to different types of matter fields, its coupling 
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to ordinary matter can be asymptotically tiny (as to satisfy constraints from gravitational 
experiments ||15||), and much stronger (as first suggested in |2"Bfl) to typical dark matter 



candidates, such as the axion. In that case, the dilaton to dark-matter coupling leads to 
an initial evolution, which is similar to the "tracking" regime [^4[] of conventional models 
of quintessence, but takes place before potential energy becomes appreciable. Later on, 
the interplay of the dark-matter dilatonic charge and of the dilaton potential leads to an 
accelerated expansion in which the relative fraction of dark energy and dark matter remains 
fixed (and of order 1), thus offering a possible explanation of the cosmic coincidence, as we 
will illustrate through explicit examples. 

The paper is organized as follows. In Sect. we present the effective string cosmol- 
ogy equations, in the small curvature -but arbitrary coupling- regime, with generic matter 
sources non-minimally coupled to the dilaton. In Sect. | IT] we discuss analytically a possi- 
bile late-time attractor characterized by a constant positive acceleration and a fixed ratio 
of dark matter and dark energy. In Sect. [IV] we provide a semi-quantitative description 
of the previous phase, during which the dilaton potential can be neglected. This phase is 
characterized by a "focusing" of the energy densities of the various components of the cos- 
mological fluid (which occurs before the epoch of matter-radiation equilibrium), and by a 
subsequent "dragging" regime in which the dilaton energy density tends to follow that of 
non-relativistic (dark) matter. We also discuss here the main phenomenological constraints 
that have to be imposed on the scenario. In Sect. [V] we consider a typical example of string 
cosmology model including radiation, baryonic and cold dark matter, and we present the 
results of explicit numerical integrations. Our conclusions are summarized in Sect. [VTl 



II. COSMOLOGICAL EQUATIONS IN THE STRING AND EINSTEIN FRAMES 



Our starting point is the string- frame, low-energy, gravi-dilaton effective action JT7J , 
to lowest order in the a' expansion, but including dilaton-dependent loop (and non- 
perturbative) corrections, encoded in a few "form factors" ip(4>), Z(<f)), a(4>), . . ., and in 



an effective dilaton potential V((j)) (see also flT5|). In formulae: 



S 



16n 



d 4 x — —F 2 + T m (4>, 



a(<j>) 



e -M)R + Z{(f>) (V0) : 
g, matter) 



Ml 



v(4>) 



(2.1) 



(conventions: metric signature: (+,—,—,—), R^a 13 = c^I\, a ^— R^ v = R atiu a ). Here 
M~ l = X s is the fundamental string- length parameter, g is the sigma- model metric minimally 
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coupled to fundamental strings, R, V are the curvature and the covariant derivative referred 
to g, and F^ u is the gauge field of some fundamental (GUT) group (a(4>) is the corresponding 
gauge coupling). We imagine having already compactified 6 dimensions and having frozen 
the corresponding moduli at the string scale. Following the basic proposal made in ||21||, we 



shall assume that the form factors i/j{4>), Z(<p), a(4>) approach a finite, physically interesting 
limit as <p — > +00 while, in the same limit, V — > 0. 

The fields appearing in the matter action T m are in general non-minimally and non- 
universally coupled to the dilaton (also because of the loop corrections [13]). Their gravita- 
tional and dilatonic "charge densities", T^ v and a, are defined as follows: 



5T 



Sg 



2 V 9 , 



8T 



9 °, 



(2.2) 



and it is important to stress that, when cf^O, the gravi-dilaton effective theory is radically 
different from a typical, Jordan-Brans-Dicke type model of scalar-tensor gravity | 25| . We 
shall give a prototype form of T m in the following section, after passing to the Einstein frame. 

The variation of Q2.1D with respect to g^ v then gives the equations: 



+ 2^" 



fJ. v Vi 
2 



e*Z -i)' 2 +<] V^V v (f) 
fa' 2 - - e^z) (V0) 2 - 2^'(V 2 0) - e^V{4>) 



/IV J 



(2.3) 



where Guv is the Einstein tensor, and a prime denotes differentiation with respect to <p. The 



variation with respect to 0, using the trace of eq. ( |2.3|) to eliminate R, leads to the equation 



fa' 2 -2e^z) (V 2 0) + 
(2i/)'V + V) + A 2 ^ 



(ZiP'-Z')+iP' fa"-3f 2 )} (V0) 



'T + a 



(2.4) 



We shall assume an isotropic, spatially flat metric background (appropriate to the present 
cosmological configuration), and a perfect fluid model of source. In the cosmic-time gauge 
we thus set 

= diag (l, -a 2 (t) <^J , = diag (p, -p 8^ , <f> = <f>(T), a = a(t), (2.5) 

and one can easily check, combining the above equations, that the matter stress tensor is 
not covariantly conserved (even in this frame), but satisfies the equation 

P+3H(p + p) = ^<p. (2.6) 
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For the purpose of this paper, and for an easier comparison with previous discussions 
of the quintessential scenario, it is however convenient to represent the dynamical evolution 
of the background in the more conventional Einstein frame, characterized by a metric g^ v 
minimally coupled to the dilaton, and defined by the conformal transformation g^ v = c\g^ v e^ . 
Here c\ parametrizes the asymptotic behaviour of ip{4>), 



c 2 = lim exp{— if) ((/))}, 



(2.7) 



and thus controls the asymptotic ratio between the string and the Planck scale, Mp = c\M 2 . 
In the Einstein frame the action (|2.1| ) becomes: 



S = -— p 1 ,l 



d X\J — g 



R- 



W,^2 2 



ilSvr / d4x ^^ F ^ + T rn(<P,c 2 i9^ i ', matter) 



(2.8) 



where we have defined 



k 2 (<p) = 3f 2 - 2e^Z, 



V = cfe^V . 



(2.9) 



For later use it is also convenient to define a canonical dilaton field by: 

,7 M P 



V2 



k((f>)d(f> 



(2.10) 



although, in solving the equations, it will be easier to work directly with the original field <fi. 

We now choose, also in the Einstein frame, the cosmic-time gauge, according to the 
rescaling 



Cl ae^ 2 , 



dt 



die*' 2 . 



a 



(2.11) 



From the (0,0) and components of eq. (|2.3| ) we obtain, respectively, the Einstein 

cosmological equations (in units such that Mp = c\M 2 = (87rG) _1 = 2) 



6H 2 = p + 
AH + QH 2 = 



-v ~ 



(2.12) 
(2.13) 



while from the dilaton equation ( |2.4| ) we get 
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k 2 (<p) (0 + 3^0) + k(4>) k'{4>) 2 + v'{4>) + l - W{(j>){p - 3p) + a] = 



(2.14) 



In the above equations H — a/a, a dot denotes differentiation with respect to the Einstein 
cosmic time, and we have used the definitions: 



h 2 ^ 2 + v(4>), n = \k\m 2 - H®. 



(2.15) 



The combination of equations fl2.12|) -( p.l4| ) leads finally to the coupled conservation equa- 
tions for the matter and dilaton energy density, respectively: 



p + 3H(p + p) - ~0 W{<t>){p - 3p) + a] = 0, 
P4> + 3H{ P<i> + P<j> ) + ^0 W{4>) (p-3p) + a} = 0. 



(2.16) 
(2.17) 



For further applications, and for a more transparent numerical integration, it is also 
convenient to parametrize the time evolution of all variables in terms of the logarithm of 
the scale factor, \ — ln(a/oj), where corresponds to the initial scalef], and to separate the 
radiation, baryonic and non-baryonic matter components of the cosmological fluid by setting 



P = Pr + Pb + Pd = Pr + Pr, 



a = a r + a b + a d = a r + a m . 



[2.18) 



The dilaton equation and the Einstein equation ( [2.12j ) can then be written, respectively, as 



2H 2 k 2 



d 



1 



dx 2 + e \2 pm + 3 pT+v )^ + 2H2kk 'y^) + 2V ' + ^ pm + = °' (2 - 19) 



H 2 



2 {d X , 



Pm+ Pr + V. 



(2.20) 



The matter evolution equation ( |2.16| ) can be split into the various components as 

0. 

dfh 

0. 
0. 



dp T ^ ^ a r d(f) 
dx 2 dx 

dpb . o 1 / , , . x d(j) 

~r + zpb - - {ip pb + cT b ) -r- 

aX ^ aX 

d Pd . q 1 , . ,d<p 
^ + 3p,--(^p, + a,)- 



(2.21) 
(2.22) 
(2.23) 



*The relation between x an d the redshift z is x = ~ l n (l + z ) + l n ( a o/ a i); where ao is the present 
value of the scale factor. 
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Finally, eq. ( 2.19|) is also equivalent to the dilaton conservation equation (|2.17 ), which 
becomes 



(2.24) 



III. ACCELERATED LATE-TIME ATTRACTORS WITH CONSTANT ft 

As a first step towards a "dilatonic" interpretation of quintessence we will now discuss 
the possibility that the above equations, together with a string-theory motivated potential 
and loop corrections, are asymptotically solved by an accelerated expansion, a > 0, with 
frozen ratio p m j ' p<j> of the order of unity. This last property, in particular, is expected to 



solve (or at least alleviate) the cosmic coincidence problem [10,11 



Under the assumption made in ^TJ that the form factors appearing in ( |2.1| ) have a finite 
limit as <fi — > +oo, and assuming the validity of an asymptotic Taylor expansion, we write: 

e -V#) = c 2 + bie -4> + {e- 2ct> ), Z(<f>) = -c\ + fee"* + 0(e-^) , 

a(ct))- 1 = a 1 + be-* + 0(e- 2<t) ) , (3.1) 

where assumed to be of the same order (typically of order 10 2 since, as already 

noted, c\ = (Mp/M s ) 2 ), and ao is to be identified with the unified gauge coupling at the GUT 



scale, i.e. ao — 1/25. Unlike the model discussed in [£6j] our model thus describes, in the 
strong coupling limit <p —> +oo, a minimally coupled, canonical scalar field <j) — V^(c2/ci)</>, 
see eq. ( |2.10|) . In the opposite limit, <fi — > — oo, the gravi-dilaton string effective action 
reduces, as usual, to an effective Brans-Dicke model with parameter u = — 1. We note that 
it is not hard to chose tp(4>) and Z(<ft) in such a way that the kinetic term of the dilaton 
keeps the correct sign at all values of <ft (see the example given in Sect. 0). 

Similarly, the assumption that V originates from non-perturbative effects, and that V — > 
as (f) — >• oo, allows us to write, quite generically: 

V{(j>) = U e-* + O(e- 2 ^). (3.2) 

Since the overall normalization of the potential Vq is non-perturbative, it should be related 
to the asymptotic value of the gauge coupling ao by an expression of the form: 



9 



with some model-dependent (one-loop) /3-function coefficient j3. For a comparison with 
earlier studies of an exponential potential P7 , |28| we also note that, when referred to the 



canonically normalized dilaton field <fi defined in ( p.lOj) , the Einstein frame potential 
asymptotically exibits an exponential behaviour V ~ exp (-\(j>/Mp), with A = ci/c 2 = y/2/k 
at — > +00. 

It is important to discuss the size of the potential needed for the viability of our scenario. 
Since the acceleration of the Universe appears to be a relatively recent phenomenon (even, 
possibly, an extremely recent one, as recently argued in p9f)> the potential V must enter the 



game very late, i.e. at an energy scale of the order of 10 3 eV. Unless we want to 

play with an unnaturally large present value of <p, this also fixes the scale of the potential 
in (|3.2|) as Vq ~ (10~ 3 eV) 4 . As far as we know, this is feature is common to all quintessence 
scenarios: the problem of an outstandingly small cosmological constant is traded for the 
introduction of another unnaturally small mass scale M*. 

In our context, we easily find that, in order to have a properly normalized potential, 
we need the constant (3 appearing in the exponent of (|3.3| ) to be somewhat smaller than 
the coefficient f3 3 of the QCD beta function (see also the discussion after eq. (|3.5|)), say 
(3 ~ 0.6/?3. Given our ignorance of the origin of the dilaton potential, this looks perfectly 
acceptable, a priori. However, this apparent resolution of the fine-tuning problem should not 
hide the fact that the potential has to be adjusted very precisely if one wants to start the 
acceleration of the Universe not earlier that at red-shift z ~ 0(1), and not later than today. 
To the best of our knowledge there is no obvious explanation, at present, of this aspect of 
the coincidence problem. 



Let us now come to the matter sector of the action ( |2.8| ). As a typical example of T m we 
take: 

r m (0, g, matter) = J d 4 x^g~ N [i$ + m N {<f))} N + ^J d A Xyf=g [e^{d^D) 2 - e vi<t>) fi 2 D 2 

(3.4) 

the first term representing baryonic matter, the second (scalar) cold dark matter, while the 
gauge term appearing explicitly in ( |2.8| ) can already represent the radiation component of 
the cosmic fluid. 

The non-observation of appreciable cosmological variations of the coupling constants 
as well as the precision tests of Newtonian gravity |3l| in the context of long-range 



dilatonic interactions, force us to assume that ordinary matter and radiation have nearly 
metric couplings to g^, i.e. that 0^, a r ~ as <p ~> °°- It is n °t hard to see how such a 
near-vanishing of dilatonic charges can be achieved starting from the actions 
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Following ref. \TM we have 



~ ~ JiZ ( ln A QCd) , — ~ S- (In a) . (3.5) 

Pb 0<P Pr OP 

Given that Aqcd ~ M s exp (— 1/ foot) (with /3 3 the coefficient of the QCD /3-function) , and 
using (|Q| ) for a, it is clear that both a h and <7 r are exponentially suppressed at large, positive 
0. This decoupling mechanism is similar in spirit to the one proposed in fl5fl , although it is 
supposed to occur at infinite bare coupling. 

In the dark matter sector, on the contrary, we shall assume more generic quantum cor- 
rections. By taking for instance the action in eq. ( |3.4|) , one has for the dilatonic charge of 
dark matter: 



0d 



C'(0) e fW (<9 M Df + r/(0) e r > w p 2 D 2 . (3.6) 



Furthermore, the equations of motion for the D field give a relation between the time- 
averaged quantities, e^^{D 2 ) = p 2 e v ^(D 2 ) (which is consistent with the interpretation of 
D as non-relativistic matter, p d = 0, as assumed in the preceding section), and relate a d and 
pd by a (generally ^-dependent) proportionality factor 

a d /p d = q{4>) = r/(0)-C'(0). (3.7) 



The lat-time behaviour we will discuss takes place if we assume that, in the strong 
coupling limit (i.e., > 1), g(0) tends to a positive constant of order unity, and that the 
dark matter component dominates over baryonic matter and radiation. Thus, the regime we 
are considering is characterized (according to eqs. ( |3.1| ), ( |3.2| )) by 

k 2 {4>) = 2cl/c\ = 2/A 2 , a = a d , p = p d , q{4>) = q = 0(1), o d = qp d . (3.8) 

It follows that the dilaton coupling to the stress tensor can be asymptotically neglected with 
respect to the coupling to the dilatonic charge, as ip' — e /ci ^ 1- The dilaton and dark 
matter conservation equations ( |2.23|) , ( |2.24| ) and the Einstein equations ( |2.12j ), ( ^.13[ ) can 



then be written, asymptotically, in the form 



q q 
p d + 3H p d - - p m (f> = 0, p^ + 6Hp k + - p m <f> = 0, (3.9) 



2H 

l = Q d + Q k + Q v , l + — = {l v -tt k , (3.10) 

6H Z 
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where we have defined 



Pd = 6H 2 £l d , p (j> = p k + pv, 

p k = 6H 2 n k = <j) 2 /\ 2 } p v = 6H 2 n v = V. (3.11) 

We now look for solutions with asymptotically frozen dark-matter over dark-energy ra- 
tio, and frozen "equation of state". From the constraint Q3.10D this is equivalent to the 
requirement that p k , pv and p d scale in the same way, i.e. 

dfogpt _ d\ogp d dlogpy _ d\ogp d 

dx dx dx dx 

The first condition and the conservation equations give 

d A = ® {Qv _n k) . ( 3.i3) 
dx q 

Expressing d(p/dx through Q k = (d(f>/dx) 2 /6^ 2 , and inserting it in the second condition 
( ft.l2j) , we obtain, respectively, 

Xq ^ ti n v -n kh , = 2^^, (3,4, 



where in the latter the asymptotic form of the potential ( |3.2| ) has been used. The last two 
equations can be solved for Qk and fly, 

giving easily 

12 + g(g + 2)A 2 q(q + 2)A 2 

= (q + 2)*\* ' W * = -12 + q (q + 2)\*> (3 - 16) 

where the last equation for = (0& — Qv)/(^k + &v) provides the dilaton's equation of 
state. 



The above asymptotic solution, first obained in H32| , and recently studied in [[33|.[|l 



generalizes the results discussed in ref. [^7] to the interacting dark matter case, and is very 
similar to the results obtained by including suitable non-minimal couplings in a Brans-Dicke 
context , or by including an effective bulk viscosity in the dark matter stress tensor 0,^] . 
Our eq. (|3.9|) corresponds indeed, formally, to a dissipative pressure II = —qp m (<f)/6H) (in 
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the notation of |p5|). See also |33| , M for a discussion of the parameter values compatible 
with such an asymptotic solution. 

Once flk and Qy are given, one can easily compute all the relevant kinematic properties 
of the asymptotic solution as a function of only two parameters, q and A = c\/c2, which 
are in principle calculable for a given string theory model. The asymptotic value of the 
acceleration, in particular, is fixed by eq. ( |3.1(J| ) as 

n H q - 1 

1 + T^ = ^T- (3-17) 



aH 2 H 2 q 

One can also easily obtain, through a simple integration, the asymptotic evolution of the 
Hubble factor and of the dominant energy density, 

#^ a -3/(2+ 9 ) ; ^^6/(2+^ (31g) 



In order to illustrate the range of parameters possibly compatible with present phe- 
nomenology, we have plotted in the {A , q} plane various curves at = Qk + Qy = const, 
and = (Qk — £V)/(^fc + ^V) = const. (Fig. 1). Note that the case discussed in Ref. 
p7| corresponds to staying on the A axis. In that case, the critical value of A below which 
VLd/VL^ — > is v^3- The addition of q makes parameter space two-dimensional, with the point 
A = replaced by the left-most curve VL^ = 1. Beyond that curve, i.e. for A 2 < 6/(2 + q) 
(as well as for all values of A if q < —2), the ratio Q^/Q^ goes to zero. However, while in the 
case of |T/J acceleration and a finite ratio Q^/Q^ are incompatible, this is perfectly possible 
in a large region of the {A , q} plane. 

In fact, it is possible to determine the region of our parameter space that survives the 
various observational constraints (Type la supernovae, CMB anisotropies, large-scale struc- 
ture . . . ). The present values of and have to lie in the range p6|j37[] 0.6 < < 0.7, 
and —1 < W(f, < —0.4, but the two allowed intervals are not uncorrelated. Assuming that we 
are already in the asymptotic regime, the allowed region lies roughly between the two curves 
fl^ = 0.6 and fl^ = 0.7 and above q = 2. Other phenomenological (but somewhat more 
model-dependent) constraints on q and A can be obtained from the recent measurements of 
the position of the third anisotropy peak in the CMB distribution |38| , which constrains the 
value of today and at last scattering, as well as the time-averaged equation of state (w^) 
||39|| . In the final part of this paper we shall present a model of dark matter that seems to 
be compatible with all the above-mentioned constraints. 
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X 

FIG. 1. The asymptotic configurations in the plane {A ,q}. The full bold curves corre- 
spond to asymptotic solutions with fixed ratios p$j Pd an d with the following values of fi^: 
1, 0.8, 0.7, 0.6, 0.5, 0.4. On the right vertical axis we have reported the corresponding q-dependent 
acceleration parameter, 'da /a 2 . The thin dashed curves correspond to fixed asymptotic values of 
the dilatonic equation of state = p^/p^, respectively —0.4, —0.7, —0.9 and —0.95. 



IV. FOCUSING AND DRAGGING WITH V = 0: AN ANALYTIC STUDY OF 

EARLY-TIME EVOLUTION 

Having discussed, in the previous section, the late-time accelerated expansion caused 
by the interplay of the dilaton potential and the dark-matter dilatonic charge, it looks 
appropriate to illustrate the earlier evolution, i.e. before the dilaton potential starts entering 
the game. In this section we shall provide a semi-quantitative, analytic analysis of this 
behaviour as it follows from the string cosmology equations ( 2.21 )— ( 2.24|) , by imposing on the 



non-perturbative normalization (|3.3|) the constraint V^ 4 <C H C<1 , where H eq is the curvature 
scale at the epoch of matter-radiation equality. In such a way the dilaton potential may 
eventually become important only at late times, in the matter-dominated era. We will show 
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that this early evolution can be roughly divided in three epochs, providing, altogether, an 
intermediate attractor that nicely connects to the accelerated behaviour described in Section 



III. 



Let us start by considering an initial, post-big bang and post-inflationary regime of 
expansion driven by the standard radiation fluid, with negligible dilatonic charge, a r = 0. 
Possible non-relativistic matter, if present, is highly subdominant with respect to the other 



components (p m <C p r ) and, consequently, the dilatonic terms in eq. (|2.24|) can be 
neglected. The conservation equations can be easily integrated to give 

p r = p ri e" 4x , p = p^e" 6 *. (4.1) 

Therefore, the dilaton (kinetic) energy density, even if initially of the same order as p r , is 
rapidly diluted like a -6 . The dilaton itself, starting from a value (fti ~ ^(1) typical of the 
moderately-strong coupling post-big bang epoch, tends to settle down to a constant value 
(as already noticed in J4T|), that can be easily estimated as follows 

»=M£V=£(£Y (4-2) 



2 \d X 12 V d Xj 



For p ri = p^ we get 



which, for k = const., leads to a solution with asymptotic value (ft = <f>\, related to the initial 
value (fti = 0(0) by the constant shift 



Atft = (ft, - & = ^ln(l + y/2) ~ | ~ -^3, (4.4) 

independently of (fti and of the initial x (the last equality holds for (fti large enough to justify 
the asymptotic relation k = a/2/ A). 

Such an initial regime is effective until the dilaton kinetic energy becomes of the same 
order as p m . At that point, some oscillations are triggered by the interference term of eq. 
O), but the dilaton energy density keeps decreasing, on the average, until it enters a 
"focusing" regime, during which it is diluted at a much slower rate (like a -2 ), so as to 
converge, at equality, towards the larger values of p m and p r . Eventually, when dark non- 
relativistic matter becomes the dominant source (pd > p r ), the dilaton energy density tends 
to follow the dark matter evolution, as if it were "dragged" by it. 
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Before turning to a quantitative analysis of these two regimes we note that the time evo- 
lution of p^, in the "tracking" quintessence, is determined by the slope of the potential. In 
the present context, instead, the focusing and dragging effects are not due to the potential, 
but they are controlled by the non-minimal coupling induced by (ip f + q) (thus implementing 
an attractor mechanism already proposed for a class of non-minimal scalar-tensor models 
of quintessence [filj). Thanks to the focusing effect, which seems to be typical of the string 
effective action (even if similar, in a sense, to the "self-adjusting" solutions of general rel- 
ativity with exponential potential [p7|), the dilaton energy density at the matter-radiation 
equality turns out to be fixed independently from its initial value, and only slightly depen- 
dent from the initial value of the dilaton, fa. For large enough values of q, however, even the 
dependence upon fa tends to disappear, because the value of the dilaton itself gets focused, 
as will be discussed in the next section. 

For a quantitative analytical study of the "focusing" and "dragging" regimes, we start 
from eqs. fl2.21|) -( pi.24 ). Lumping together baryonic and dark matter, neglecting V, and 



assuming, according to (|3.7|), a = a rn = q(fa) p m , those equations can be easily recast in the 
form: 

p; 1 ^ + 4 = 0, (4.5) 
Pm 1 ^+[3TV3e(p,/p)^]=0, (4.6) 
+ 6pt±V3e Pm (p (t> /p) 1 / 2 = 0, (4.7) 



dx 

where we have introduced the important parameter: 

e(0) = k(<f>) ' (4 - 8) 

and the sign ambiguity comes from solving eq. (|4.2j ) for d(f>/dx in terms of p^. The focusing 
solution is then characterized by the relation: 

„ _. n 2 {fa)pl 



(4.9) 



i.e. = n 2 (fa) Q 2 m , which holds under the assumption that both e and n are slowly varying. 
Indeed, we can establish the connection between these two quantities by inserting the ansatz 
( PI) into 0. This gives: 
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- 6T V3i = 2„-^ + 2p-'^-p-'* (4.10) 
n dx dx dx 

where on the right-hand side the logarithmic derivative of ( f4.9|) has been taken. By using 
( [4.6|) one finally has 

p- 1 ^ T ^3 e \ n - 1 + 2tt m n] = 2rr 1 -^ ~ 0. (4.11) 
dx 1 dx 

We can now discuss a few cases of interest. During the radiation-dominated phase, and 
after the kinetic energy of the dilaton is quickly red-shifted away, we can neglect the term 
with Q m in eq. fl4.11|) , we set dp/dx = —4p, and obtain: 



v^e d(f) 3p m e 3/4e 2 

n ~ , — ~ : — , p& ~ — - — . 4.12) 

4 ' d X 2k p ' H4> 16 p K J 

We refer to this behaviour as "focusing" since it implies that p m lies, modulo a factor 
(16/3)e -2 , at the geometric mean between p ~ p r and p^. Hence, as we approach radiation- 
matter equality, p^ is effectively focused towards the same common value of the other two 
components (see eq. ( |4 . 1 5| ) below). Note that, for a positive e, this happens thanks to a 
negative d(j)/dx- 

In the matter-dominated regime it is no longer safe to neglect the term in Q m in eq. 
( |4.11| ), unless eC 1. In that case, the solution is 



n -H- p ~tT' n ^ 3 ^(- p ~ty 2 ■ <4i3) 

During matter domination, using dp/dx = ~ 3p, one gets 

e d(f> 2e p m e 2 

p **—- (414) 

In other words, the focusing regime has been turned into a dragging one: the dilaton energy 
is dragged along by the (dark) matter energy and keeps a (small) constant ratio to it. 
Incidentally, at the epoch of exact matter-radiation equality, using dp/dx — ~ 3.5/0, we 
easily get (still at small e): 

^ ~ — , (4.15) 

Pe q 49 V ^ 

which is always smaller than 6% for e < 1. 
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In order to understand what happens at larger values of e it is useful to find the reason 
why, for small e, p^j p m stays constant. This comes about because the corrections to the a~ 3 
and a -6 laws for p m and p^, due to the non-vanishing e, push the two towards each other. 
It is easy to check that, precisely if p^/{p m + p<j,) = e 2 /3, both energies scale like a~^ +e2 \ 
We note, incidentally, that the above ratio of energies nicely fits with the value given in 
( 4.14 ) when e < 1. If e < 1, the decrease of p^ is still slower than the a -4 of p r , which 



justifies neglecting the latter. However, if e > 1, this is no longer the case and we have a 
third kind of behaviour, which can be called "total dragging" . In that shown by a 

simple analysis, all three components of p scale like radiation, with the following sharing of 
the "energy budget" (remember that we are always at Q = 1): 

n, = £ = ^, n, = ^i. (4.i6, 

In the next section we will see how numerical integration confirms in full detail the analytic 
behaviour we have discussed. We end this Section by discussing some constraints on our 
parameters. 

As already mentioned, we assume the ordinary components of matter (radiation and 
baryons) to have a nearly metric coupling to g^ u (see discussion after eq. (|3.4j) ). To be 
more specific, let us define the ratios between dilatonic charges and energy densities in a 
way similar to that used for cold dark matter in eq. (|3.7|), i.e. 

q r (4>) = a r /p r , q b (4>) = a b /p b . (4.17) 

Since it is precisely the ratio (ip' + q r ,b)/k, which controls both the effective coupling of 
the dilaton to ordinary macroscopic matter, as well as a possible time-dependence of the 
fundamental constants p2| , ^3| , we shall assume that both q b and q r are at most of order 



ip', in agreement with the discussion after eq. ( |3.5|) . We then find that there are neither 
appreciable violations of the equivalence principle in the context of macroscopic gravitational 
interactions, nor significant contributions to the time- variation of the fundamental constants, 
both effects being controlled by ip'/k for q r ^ —> 0. In the strong coupling regime we have 
if}'/k ~ e. J{c\C2). For a non-negative <pi, and c\,c\ of order 10 2 , there is no appreciable 
deviation from the standard cosmological scenario down to the epoch of matter-radiation 
equality, so that one easily satisfies the early- Universe constraints on dark energy, as reported 
for instance in |44[ . 



The dilaton charge of dark matter is not restricted by the experimental tests of long- 
range gravitational interactions: this is the reason why we can play with it in order to 
produce an acceleration. Still, from the above discussion on the early phases of the universe, 
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it is clear that high values of the dark-matter parameter e may result in dangerously high 
values for fl^, and thus in radical deviations from the standard cosmological scenario. Until 
radiation-matter equality the situation is relatively harmless: we can easily estimate the 
dilaton energy density at the equality and at the nucleosynthesis scale, Hn ~ 10 10 H i 



eq • 



using the fact that the dilaton, during the focusing regime, is not significantly shifted away 
from the value fa + A<f>, fixed by eq. ( |4.4| ). Because of the focusing behaviour we find 
f^(nucl) ~ 10~ 10 f^(eq), and therefore the most stringent bound comes at equality, where, 



thanks to ( 4.15|) , it is comfortably satisfied for e < 1. 

During the dragging phase, however, we must certainly impose e < 1, otherwise, the 
phenomenon of "total dragging" takes place. This would represent a dramatic deviation 
from the standard cosmological scenario, since all the components p^, p r , pd (except baryonic 
matter) would redshift in the same way (a~ 4 ) from equality until the potential starts to be 
felt. Even if e < 1, but not sufficiently small, the unusual scaling p m oc a~ 3 ~ e2 tends to 
change the global temporal picture between now and the epoch of matter-radiation equality 



and, from eqs. ( (4.14| ), values of 0,$ ~ e 2 /3 (while in agreement with possible constraints at 



last scattering f44]| ) can be dangerously high. In our context, a bound fi^(drag) < 0.1, i.e. 
e 2 (drag) < 0.3, appears to be necessary in order to agree with the observed CMB spectrum 
and with the standard scenario of structure formation. 

On the other hand, due to the smallness of ip' ~ e~^/cf in the dragging regime, an upper 
bound on e effectively turns into a bound on the value of q/k, i.e. on the dilatonic charge 
of the dark matter component. The above constraints thus translate into a bound for the 
combination Xq: 

Ag(0dra g ) < 0.8, (4.18) 

where we used the already mentioned asymptotic relation A = y/2/k. It is clear that a 
constant q cannot satisfy the above bound and, at the same time, provide the present 
acceleration of the Universe by means of the mechanism described in Section |l| (see also 
Fig. 1), that requires q\ > 4. 

A time- (or, better, 0- ) dependent q, however, is allowed. For this reason we have 
to consider cold dark matter models like the one of eq. ( |3.4|) , whose dilatonic charge ( |3.7| ) 
switches on at large enough values of the dilaton. The transition to large values of is 
rapidly activated as the potential comes into play, pv ~ P4>- At that point, the dilaton 
energy density stops decreasing and freezes at a constant value, necessarily crossing, at some 
later moment, the matter energy density, p^ ~ p<i- From then on, the dilaton starts rolling 
towards +oo, triggering the effect of the dilatonic charge, which rapidly freezes the ratio 
p(l>/pm and (for suitable values of q) leads to the accelerated asymptotic regime described 
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by eqs. ( |3.16|) -( j3.17|) . Explicit numerical examples of such a behaviour will be discussed in 
Section M. 

For a realistic picture, in which the positive acceleration regime starts around the present 
epoch (and not much earlier) and the standard scenario of structure formation is implemented 
successfully, we have to require that the contribution of the dilatonic charge (as well as the 
effect of the dilaton potential) come into play only at a late enough epoch. The importance 
of this constraint was already discussed in the context of other scalar-tensor models of 
quintessence [28| where, for instance, the non-minimal coupling of the scalar field to the trace 



of the dark matter stress tensor was assumed to be ^-dependent, to interpolate between a 
small and a large mixing regime. 

V. NUMERICAL EXAMPLES 

Finally, after the analytic discussion of the previous section, it seems appropriate to 
illustrate the "run-away" dilaton scenario with some numerical example, both in order to 
confirm the validity of some approximations made in deriving the analytic results, and in 
order to see how the various regimes we discussed can be put together. To this aim, we shall 
numerically integrate eqs. ( |2.19|) -( p.23| ), using eq. (|2.20| ) as a constraint on the set of initial 



data, and assuming an explicit model for the dilatonic charges and the dilaton potential. 



Also, following the "induced-gravity" ideas [21], we shall specialize the loop form-factors 



according to eq. (|3.1| ), using the "minimal" choice 

e -m = e -* + c 2 1} Z(<f)) = e~*-cl. (5.1) 

First of all, for a clear illustration of the "focusing" and "dragging" regimes, let us put 
V = 0, a r = = crft, and ad = qpd, with q = const. By choosing, in particular, c\ = 100, 
c 2 . = 30, we have integrated eqs. ( |2.19| )-( p.23| ) for three different values of the charge, q = 0, 
q = 0.01, and q = 0.1, starting from the initial scale Hi = 10 40 H cq , 

,^ 1/2 = (fhm\ = in - 20 

V H 




/ TT ,1/2 

' _//^M_ 10 -m (52) 




and using p^i = p ri , (pi = —2 as initial conditions. It should be noted that such initial 
conditions are generic, in the sense that different initial values of p^ and may change the 
fixed value reached by <fi during the focusing phase, but do not affect in a significant way the 
subsequent evolution, as will be discussed at the end of this section. 
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FIG. 2. Time evolution of for q = (dash-dotted curve), q = 0.01 (dashed curve) and 
q = 0.1 (dotted curve). The initial scale is a, = 10 _20 a e q, and the epoch of matter-radiation 
equality corresponds to x — 46. Left panel: the dilaton energy density is compared with the 
radiation (thin solid curve) and matter (bold solid curve) energy density Right panel: the dilaton 
energy density (in critical units) is compared with the analytical estimates ( \4.1ty , (4.13), ( 4.14\) for 
the focusing and dragging phases. 



The results of this first numerical integration are illustrated in Fig. 2. The left panel 
clearly displays the initial regime of fast dilaton dilution (p^ ~ a -6 ), the subsequent focusing 
regime (p^ ~ a -2 , see eq. ( |4.12| )) triggered (after some oscillations) soon after p^ falls below 
p m , and the final dragging regime (p^ ~ p m , see eq. Q4.14D ) in the epoch of matter domination 
(the epoch of matter-radiation equality corresponds to x — 46). Note that the constant 
values of q have been chosen small enough to avoid the phenomenon of "total dragging", 
see Section |V|. Note also that, in this example, p m always coincides with pd- In the right 
panel the evolution of f^, obtained through the numerical integration, is compared with the 
analytic estimates ( |4.12| ), ( f4.13| ), ( |4.14| ), for the three different values of q. In all cases, Q^, 
grows like a 2 during the focusing regime (in the radiation era), while the final stabilization 

= const, after the epoch of matter-radiation equality (% > 46), clearly illustrates the 
effect of the dragging phase during which p^ and p m evolve in time with the same behaviour. 

For a realistic model of quintessence, however, a constant dilatonic charge cannot drive 
the Universe towards an asymptotic accelerated regime and, simultaneously, satisfy all the 
required phenomeno logical constraints during the earlier epochs (as discussed in the previ- 



21 



ous sections). By keeping o&, a r ~ at large coupling (see eq. (|375| ) and the discussion 
thereafter), we shall thus consider the explicit model of scalar dark matter fl3.4|) , with the 
following simple loop form-factors 

e -CW = i + e <?o</y c 2 ; e m = const (5.3) 

(note that, by a field redefinition, one of the two loop factors can always be taken to be 
trivial: what really matters is the ratio e^/e v ). Using ( |3.7| ) we immediately get 

«W=^ = *?T^' (5 ' 4) 

which is exponentially suppressed in the perturbative regime, and approaches q = q$ sX large 
coupling (for q Q > 1 it is thus compatible with an asymptotically accelerated cosmological 
configuration, see Fig. 1). For our numerical example we shall choose go = 2.5 and c 2 = 150, 
but the behaviour of the solution is rather stable, at late times, against large variations of 
the latter parameter (see the discussion at the end of this section). 

In addition, we have to specify the form of the dilaton potential. In agreement with 
its non-perturbative origin, and with the assumtion of exponential suppression at strong 
coupling (see Section |TV| ), the simplest choice is a difference of terms of the type e~Pl a ^ . 
We shall thus consider the bell-like potential (in units Mp = 2) 

V{<j>) = m 2 [exp (-e-Vft) - exp (-e^/ft)! , < (3 2 < fa, (5.5) 



which leads, asymptotically, to the large-0 behaviour of eq. ( |3.2|) . The mass scale my, related 
to the mass M* of eq. (|3.3| ), will be fixed at my = 10~ 3 H cq , together with (3\ = 10, /3 2 = 5, 
for a realistic scenario that starts accelerating at a phenomenologically acceptable epoch. 

With all the parameters fixed, we have numerically integrated the evolution equations 
( 2.19 )— ( 2.23 ), for our model of charge ( |5.4j) and potential (|5.5|) , using the same initial con- 



ditions as in the previous example, but separating the dark and baryonic components inside 
p m . In particular, we have set, initially, pdi = 10~ 20 p ri , p bi = 7 x 10~ 21 p„. 



The resulting late-time evolution of the various energy densities is shown in the left 
panel of Fig. 3. Dark matter and baryonic energy densities evolve in the same way, until 
the potential comes into play, starting at a scale around x — 49. The potential first tends to 
stabilize to a constant but then (thanks to the contribution of q) the system eventually 
evolves towards a final regime in which p^ and pd are closely tied up, and their asymptotic 
evolution departs from the trajectory of the standard, decelerated scenario (in particular, 
they both scale, asymptotically as a~ 6 ^ 2+qo \ see eq. (|3.18| ). It is amusing to conjecture that 
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FIG. 3. Left panel: Late-time evolution of the dark matter (solid curve), barionic matter 
(dashed curve), radiation (dotted curve) and the dilaton (dash-dotted curve) energy densities, 
for the string cosmology model specified by eqs. ( 5.4), ( |5.5| ). The upper horizontal axis gives the 
log 10 of the redshift parameter. Right panel: for the same model, the late-time evolution of q 
(hne-dashed curve), (dash-dotted curve), (solid curve) and of the acceleration parameter 
'da /a 2 (dashed curve). 



the different time-dependence of pb and could be responsible for the present small ratio 
Pb/pd- 

In the right panel we have plotted the time evolution of the dilatonic charge q, of the 
energy density fl^, of the equation of state w^, and of the acceleration parameter d/aH 2 . 
When the potential energy becomes important, all the above quantities rapidly approach 
their asymptotic values given in eqs. (|3.16|) -( 3.17 ). Note that, with our choice of parameters, 



we have go = 2.5 and A = C1/C2 = ylO/3, corresponding to an asymptotic value ~ 0.733, 
slightly exceeding the best fit value suggested by present observations [^,0. It is important 
to stress, however, that the asymptotic attractor may be preceded by a (short) oscillating 
regime, which, as illustrated in the right panel of Fig. 3, can easily allow for values of the 
cosmological parameters different from the asymptotic ones to be compatible with present 
observations. Note also that, when switching from the focusing to the dragging phases, the 
dilaton starts to move back towards decreasing values of q, as will be illustrated also by a 
subsequent numerical integration. This may slow down the evolution of p^ with respect to p m 
during the dragging, as shown for instance in the left panel of Fig. 3. Because of this effect, 
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FIG. 4. Time evolution of the dilaton Geld, for different initial conditions (pi = —4, —2, 0, 2. A7i 
tie other parameters are the same as in the example of Fig. 3. After the plateau associated with 
the focusing regime, and for a strong enough dilatonic charge, the solutions tend to converge to a 
common value of 4>. The subsequent running to +oo, driven by the potential, is thus completely 
independent from the initial value. 

however, the dilaton can easily satisfy, during the dragging phase, the phenomenological 
bounds discussed in the previous sections. This does not require fine tuning, the validity if 
the bounds being guaranteed for a large basin of initial conditions by a convergent behaviour 
of the solutions during dragging. 

During the focusing phase, in fact, the dilaton is practically frozen, as can be argued 
from eq. Q4.12| ), and its effective constant value, as determined by eq. (|4.4|), depends on 
(pi. However, if such a value is high enough, the presence of the dilatonic charge may 
become important, and may contribute to the focalization towards the epoch of matter- 
radiation equality, as already anticipated. This is illustrated in Fig. 4, which shows the time 
evolution of the dilaton obtained by numerically integrating the same model as in Fig. 3, 
for different initial values <pi = —4, —2, 0, 2. Although we start with different dilaton values 
at the plateau associated with the focusing regime, all the solutions tend to converge as we 
enter the dragging regime, so as to make the subsequent (potential-dominated) evolution 
insensitive to the initial value of the dilaton^. 



^The preceding evolution, of course, is not sensitive either, since during focusing the order of 
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FIG. 5. Time evolution of q{4>), from eq. (5.4), for three different values of the parameter c. 
All the other parameters are the same as in the example of Fig. 3. During the dragging phase the 
value of q converges to the regime q <C 1 . 



This new focusing effect, which is very different from the one of the energy densities 
during the radiation-dominated phase, can also be understood analytically by writing the 
solution of eq. ( f4.14|) as: 



X Xcq 



2e(0) 



(5.6) 



Since k is almost constant, a variation 5<f) cq on the initial value of <fi changes the solution 
4>{x) by an amount 8(f>(x) = [e(0)/e(0 eq )]50 e q, which rapidly decreases (with q(<f>)) during 
the dragging phase. This is why the solution has become independent of the initial value of 
by the time the potential becomes an important component. 

For the same reason, the model is only weakly affected by variations of the parameter c in 
eq. (|5.4f) , which roughly gives the transition scale between small and large dilatonic charges: 
4> s = (2/g ) logc. Indeed, because of the above mechanism the dilaton is pushed back during 
the dragging phase the dilaton is pushed back, with a velocity as high as needed to reach, 
in any case, the safe zone q <C 1. This effect is illustrated in Fig. 5, where we have plotted 
the time evolution of q(4>), for the same model as Fig. 3, and for three different values of c. 

It should be noted, in conclusion, that the above class of models depends in crucial way 



magnitude of VL^ is given by $7^ as i n 
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on three important parameters: my, q and the ratio A = C\jci- The first one controls the 
transition time between the epoch of standard cosmological evolution and the final acceler- 
ated regime (as can be easily checked, for instance, by repeating the numerical integration 
of Fig. 3 with different values of my). The other two parameters control the asymptotic 
properties of the model (acceleration, equation of state, ...), as discussed in Section |TTTj. 
Future precision data, both from supernovae observations and from measurements of the 
CMB anisotropy, could give us a good determination of these parameters, thus providing 
important indirect information on the parameters of the string effective action in the strong 
coupling regime. 

VI. CONCLUSION 

Let us conclude by summarizing the main points of our work. We have argued that a 
run-away dilaton can provide an interesting model of quintessence under a well-defined set 
of assumptions that we list hereafter: 

• The limit of supertring theory, as its bare four- dimensional effective coupling goes 
to infinity (so called induced gravity/gauge or compositeness limit), should exist and 
make sense phenomenologically, i.e. should yield reasonable values for the unified gauge 
coupling at the GUT scale and for the ratio Mp/Mqut-, thanks to the large number of 
degrees of freedom at Mqut] 

• In the visible-matter sector, the couplings to the dilaton, either direct or through the 
trace of the energy-momentum tensor (i.e. via a conformally rescaled metric), should 
vanish in the — > +00 limit; 

• In the dark matter sector, there should be a surviving coupling to the dilaton (and 
thus violations of the strong and/or weak equivalence principles) even in the (ft — > +00 
limit; 

• The dilaton potential should be non-perturbative, go to zero asymptotically, and have 
an absolute scale not too far from the present energy density. 

Under these circumstances, it is natural for the dilaton energy in critical units, Q^, to 
be: i) subdominant during radiation domination; ii) a (small) fraction of the total energy 
at matter-radiation equality; iii) a (small) fraction of fl m during the earlier epoch of matter 
domination; iv) a fraction of dark-matter energy since a red-shift 0(1). This very last phase 
is characterized by an accelerated expansion. 
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In other words, this framework seems to be naturally consistent with present astrophysical 
observations and with known cosmological constraints. From a theoretical point of view the 



model appears to combine nicely previous ideas |15j on dilaton stabilization and decoupling 
with those recently advocated (e.g. in fl28l ) so as to have acceleration while keeping the ratio 
fi(dark energy) /Q( dark matter) constant. 

It must be stressed, of course, that the analysis presented in this paper is still prelimi- 
nary, and that various problems are still open. In particular, a precise computation of the 
CMB anisotropy spectrum, and of the spacing of acoustic peaks expected in this context, 
as well as a comparison with currently available measurements , could provide significant 
bounds on the parameters of the string cosmology models we have discussed. Such an in- 
vestigation is postponed to future work. Nevertheless, we believe that the results of this 
paper are encouraging, as they suggest that the dilaton, which can hardly play the role of 
the inflaton in the standard inflationary scenario , could play instead a successful role as 
the quintessential field in post-inflationary, late-time cosmology. 
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